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NUMERICAL SIMULATION OF TURBULENT JET NOISE - PART 1 

By V--'; 

Ralph W. Metcalfe and Steven A. Orszag 

Flow Research , Inc. 

Kent, Washington 


Summar y ■ 

This report describes the numerical simulation of a low-Mach 
number turbulent jet using spectral methods with (32) = 32,768 

Fourier modes. The purpose is to gain insight into the generation 
of aerodynamic noise by an examination of certain detailed flow 
characteristics, such as quadrupole moments, which cannot yet be 
measured in the laboratory. The mean flow quantities are set in 
accordance with the experimental data of Liu & Maestrello^ 
and the incompressible Navier-Stokes equations are then solved 
numerically. Due to computer speed and memory constraints, iso- 
lated downstream sections of th e jet are modelled separately with 
the mean flow characteristics held constant in time. The flows 
are allowed to evolve until the fluctuating velocity components 
reach a statistically steady state. Cross section contour plots 
of the velocity components and the quadrupole moments at three 
different downstream positions are presented. 



1. Introduction 


Noise production by aerodynamic jets is a basic fluid 
mechanical problem with important practical consequences in 
the design of aircraft engines. In this report, numerical 
simulations of an idealized jet flow are discussed with the 
the goal of obtaining better understanding of the mechanism 
of noise production. 

Laboratory experiments, such as those reported by Liu & 
Maestrello^, can provide relatively complete information on 
simple flow quantities, like the mean velocity field, but are 
not yet capable of resolving the quadrupole moments of the 
flow which govern noise production. On the other hand, 
numerical simulation permits detailed examination of the flow 
at any point in time, and the calculation of arbitrary flow 
parameters, including quadrupole moments, at all space points 
However, numerical simulations also have their drawbacks, 
especially limited resolution which restricts the Reynolds 
number of the flows that can be accurately simulated. It 
is preferable to view laboratory and computer experiments 
as complementary, not mutually exclusive; in the present 
work, laboratory data are used to fix certain gross features 
of the jet flow, which | are then examined in more detail by 
numerical simulation. 

In Sec. 2, the flow model and the equations of motion 
are discussed. In Sec. 3, the numerical method is presented , 
while in Sec. 4 the numerical simulations performed to date 
are reviewed. Finally, in Sec. 5, a discussion is given of 
how the present results will be used in the study of noise 
production. 
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2. Flow Model 

As a first step in the numerical simulation of jets/ one 
can restrict attention to low Mach numbers., so that the Navier- 
Stokes equations for incompressible flow govern the evolution 
of the rotational components of the Velocity. At low Mach number 

the acoustic (compressive) part of the velocity couples weakly 

. " : I 2 

with the rotational components , so it is possible to determine 
the weakly excited acoustic field after having solved the dynamic 
equatibns for the fluid motion. The source terms in the wave 
equation for the acoustic component are just the quadrupole 
moments ; V-/- 

f)V. I)v. 

“ S3«\ (1 > 

; 3 :. . i 


where summation over repeated indices is implied and v is the 
three dimensional velocity field. Lighthil I 's theory implies 
that the radiated sound field (perturbation density) is propor- 


tional to 



_ 1 __ 

. 2 
4nc 


Q(y it 




•> 


x-y 


( 2 ) 


where c is the sound speed and Q is calculated in the incom- 
pressible limit (with an error of order Mach number squared) , 

The Navi or- Stokes equations for incompressible flow a re 




+ v(x,t) • Vv(x;t) - -Vp(x,t) + 


vV^v ( x, t ) 


( 3 ) 


V* v (x, t) 0 

where p is the pressure and v is the kinematic viscosity. 

The region of importance for noise production in a jet extends 
from the nozule downstream to at least: 20 jet diameters, while 
the lateral flow includes the jet together with, the surrounding 
potential flow. Quite a few difficulties arise, in the attempt 
to simulate this flow numerically. Among the most important are: 


1. Limited spatial resolution because of limited computer 

power. Even using a CDC 7600, it is hardly feasible to consider 
. . 5 

calculations with more than about 10 degrees of freedom to re- 
solve the velocity field. With a three dimensional flow field, 
this limit implies that no more than about 50 degrees of freedom 
can be used to resolve each flow direction- This is not suffi- 
cient to resolve twenty jet diameters downstream. 

2. Boundary conditions for a realistic flow field are dif- 

ficult, though not impossible, to impose. Some of the problems 
concern: a) the upstream boundary, whore the inflow is not known 

in all detail; b) the downstream boundary, where outflow bound- 
ary conditions may cause numerical instabilities; and c) the 
lateral boundaries, where fluid must bo entrained into the jot 

at the proper rate in order to maintain the proper jot growth. 

Because of these difficulties , the problem was further re- 
stricted as will now be described. As discussed above, labora- 
tory experiments have dete rmined the large-scale mean- flow 
characteristics, so that the goal of numerical oxpori men I at ion 
on jet noise should be the resolution of the quadrupo ie-momon t 
distribution. Since the quadrupolc source terms are dominated 
by small-scale motions and since the small scales of turbulent 

y v '3 : 

flows largely decouple from large scales , it follows that the 
distribution may be studied by isolating a local section of 
the jet. This suggests slicing the jet into several sections 
in the axial direction, and then, within each section, imposing 
the experimentally determined mean-flow characteristics on the 
flow. The mean velocity profile within each sectiou is held 
fixed in time during the numerical calculations, and the turbulent 
fluctuations are allowed to approach a statistically steady state. 
After evolving the flow in this way, the quadrupole moments are 
obtained and can be used in the Lighthill integral to determine 
the radiated sound field. 
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At the present time, the flow is also simplified by neglect- 
ing the entrainment of fluid at the lateral boundaries. The flow 
is assumed to take place in a rectangular box with periodic bound- 
ary conditions applied in the downstream direction (:•;,) and no- 
stress boundciry conditions applied in the' lateral directions (x 2 
and x^) . 

The velocity field is expressed as 

V (X, t) = V ( X ) + v' (x, t) 

where v is the mean velocity and v' is the fluctuating velocity 
The mean velocity v is held constant in time and is of the form 

V- ^ ^ = u (r)x 

where is a unit vector in the -di rection , r is the dis- 

tance from the jet axis { in the center of the computational 
domain), and U (r) is the experiment ally observed profile. At 
the present level of approximation, there is no mean radial velo- 
city since the downstream variation of mean velocity is so If- con- 
sistently neglected (conservation of mass V* v - 0 implies that, 
axial variation of v induces moan radial, velocity) . The proce- 
dure for handling the fluctuating velocity is discussed in Sects. 

3 and . 4. 


3, Num er ical Methods 

The Navier-Stokes equations (3), (4) can also be written in 

rotation form. 


3 v ( x , t ) *-,'>■ , * > , > . . > ' > L . , , , 2 > . > , , 

ij t~ _ “ v ? t) (>G U -Vu (x, t) +v»V v ( x , 1 ) 


(5) 


> , . •> > > y r 

where v(x,t) is the velocity vector, m (x, t) - V v -v(x,t) is the 
vorticity, 7i(x,t) = p(x,t) -t 1/2 | v ( x , t ) | ^ is the pressure head, 
p ( x , t ) is the pressure and v is the kinematic viscosity. The 

rotation form (5) is advantageous when using pseudospectral 

approximations since in this form the equations semiconserve 
4 

energy . This pointwise energy conservation prevents aliasing 
errors from inducing nonlinear instability, although it has 
been shown 5 that the appearance of nonlinear instability due 
to aliasing errors can indicate insufficient resolution. 

Given the structure of the problem as described in Section 
2, i L is convenient to use the following Fourier series roprosen- 
tat Lon for the velocity field v (x , t ) - ( v^ (x , t) , v 2 ( x , t ) , 
v_ (x, t) 


k 2 k 3 


V ( X , t ) 
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)' y u (k , t) exp f 2iri k , X- /L . J 


Ikil-Ki k 2 =0 -k-- : 0 


{ cos 

i ‘^2 X 2 // ^2 C ° S V ^ ° X 

3 /L 3 " - 

:< l sin 

ji k 2 x 2 / h 2 co s it k ^ x 

3 /L 3 “ “ 

1 cos 

"I k_x„ / L „ s i n i! k 1 x 

3 /I .3 a - 

indices k 

_ are intecjral 

and L . 
1 

direction. 

Choosing the 

summation 


m 


as K ' = 16, K = K = 32, permits the represent at- i on of each 

X t J 

velocity component by (32)' = 32768 Fourier modes so that the 
entire velocity -field involves about 1,0 independent degrees 


of freedom. The boundary conditions are periodic at x ^ = 0 , 


and are stress- free at 


x 2 = 0, 


and x. 


0, 


V 


3v, 

w. 

d 

' jv 1 

:ix. 


3v 3 

3^2 = V 2 = ° ° n *2 = " 


i v 2 

;jx„ 


v 3 = 0 on x.^ •• 0 , 


I,. 


I. 


V (x.^ + mL^ ,x 2 f X-j) - v ( x j , x,„ 1 X 3 ) 


r» integral (7) 


Equations for the Fourier components u(k) are derived 

4 

by Galerkin's procedure . After applying the incom- 
pressibility condition (4). to eliminate the pressure, the equa- 


tions for u(k) are found to bo (in the special case I, 
L 2 = L 3 = ’ ,l) 


2 >1 , 


dCi (k,t) ? 

- + vk ' u (k , t) = - k .. (6 -k k , 
a 3 try a y 


> > 


2 


>)t 


) ). 0 (p, t) i» (q , t) 

> * ,> V 2 > 

p-t q-k 

| p - I < K . 

1 1 J. 1 


( 8 ) 


q. 


where 


u ri (k, t) = 


■ i 
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4 U l (k l' 


k 3 | ,t) 


~sgn (k 2 ) u 2 (k^, [k_ 2 ! , | k 3 | , t) 


^ 4i sgn <k 3 U 3 {k ] ' i k 2 1 ' l k 3 ! i 1;) 


<K . 
1 


U = 1 
a ~ 2 
- 3 


(9) 


The system (8) can be solved numerically by using fast Fourier, 
transforms to evaluate the convolution sums, leap-frog time 
differencing on the nonlinear terms, and Crank-Nicolson (im- 
plicit) time differencing on the viscous terms. This method 
of solution requires the use of 18 Fourier transforms on (32)' 
points per time step. 


With the representation of the velocity field requiring 

5 . 'i ' "1 

about 10 computer storage locations , the current code requires 
the buffering of data between small core memory, large core mem- 
ory and disks. Small core memory on the CDC 7600 has 65K words 
and large core has 250K words. The code has been sped up by 
about a factor of 2 over the original Fortran version by pro- 
gramming many of the critical internal loops in assembly language, 
and now requires about 6 seconds of compute.: time per time step. 
The runs discussed later required from 300 to about 1000 time 
steps to achieve a statistically steady state. 


4. Description of the Simulations 

j The initial mean velocity profile in the numerical experi- 
ments was chosen to correspond to the experimentally measured 
mean velocity profile of an axisymmetric subsonic cold jet ' . 
The experimental apparatus consisted of a long pipe with a nozzle 

diameter of d^ = .0625 meters. The mean o.xit voloci ty was 
213m/ sec corresponding to a Mach number of 0.62. The cylin- 
drical coordinate system; (X^r) is used, where X 1 is th^ distance 
(expressed in units of the jet diameter) along the jet axis from 
the virtual origin of the jet and r corresponds to the radial 
distance from the jet axis. The jot spreads downs l ream, its 
boundary making an angle of 


0 .= 11 ° 
e 


with the X.. axis. In 


this flow, the virtual origin of the jet corresponds to a dis- 
tance 


Xq = ^-2.57 upstream of the mouth of the nozzle, located 


at 


X 1 - V 


The mean velocity profiles have been characterized in two 
self-similar regions: a mixing region, just downstream of the 

jet nozzle, corresponding to 0 < D < 5, ana a cully LurbulenL 

region for D > 5, where D = - Xq. In particular, for 

o < d < 5, 


U(X, , r) 

u 


x ,r) x 1, r < .5 - . 102 ( 1+D) 

( X x ) \ exp (-61.53 (ii + . 1 0 2 ) 2 ) , r > . 5- . 1 0 2 ( 1 +D ) 


(10) 


where U c (X^) is the mean centerline jet velocity , and 


_ r-._5 
y 2 1+D 


(ID 


and for D > 5 , 


U ( X x , r ) 


U C (X 1> 


~ — = exp (-1 . 44 (|i +. 69) ) 

if \ J 


( 12 ) 
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where 


y = 10. 3{~-}- .69 
x i 


( 13 ) 


Superimposed on the given mean velocity field is a turbulent 
fluctuating component v'(x, t). To ensure that v'(x,t) satis- 
fies the incompressibility constraint (4), we express v' (x,t) 
in the form 


} v* (x,t) = V x A ' (x) 

where the vector potential A' (x) is of the form 


(14) 


A' (x) = (I (r) ] 1/2 B (x ) 




Here, B (x) is chosen as a realization of a homogeneous, iso^ 
tropic random field with specified isotropic energy spectrum 
E(k), while I(r), the intensity function, is chosen in accord- 
ance with the experimental data on the root-mean-square turbulent 
intensity 1 . E(k) is of the form 


E (k) = Ak^exp (-Bk 2 ) . 

I (r) has been approximated at various downstream positions 
e.g. at D = 4 


where s = 


R(s) = [I (r) ] 
n . Also, we use 


2 

1/2 f •16-.ilexp(-5s ) 0<s‘ 1 


.159 exp (-1. 5 (s-1) ) s> 1 


R (s) = 


R(s) » 


.185 - .055exp(-8s ) 0<s< . 55 


. 183exp(-. 7(s-. 65) ) 


s> . 65 


. 225 - .045exp(-16s ) 0<s- . 4 


for D=8 


for D= 12 


. 2215exp ( - . 5 (s- . 4) ) 


s > . 4 


(16) 


D, 


(17) 


( 18 ) 


Three principal runs are reported corresponding to loca- 

.. V ■ .. •' : • ! i • : ! - 

tions D - 4, 8, and 12 jet diameters downstream of the jet 

exit. Each computational slab extends over about 2-3 jot 
diameters in the downstream direction (the variation being due 
to the need for increased size of the computational box further 
downstream). The initial mean flow profiles are shown in Figs. 
1 , 2 . \ 1 


Finally, it is important to discuss the Reynolds numbers 
of the simulations. In all three simulations; v = . 006m /s 
so that the Reynolds number based on jet diameter and centerline 
velocity is about R ~ 22,000. On the other., hand, the laboratory 

-5 2 

jet is run m air with the much lower viscosity v - 1.5x10 ~ m /s 
giving a Reynolds number roughly 40 times larger. The important 
question concerns how these two flows, numerical and laboratory, 
can possibly relate to each other considering the huge discrep- 
ancy in Reynolds number. However, this discrepancy in Reynolds 
numbers appears to be important only at the very smallest scales 


and not at those scales of motion responsible for sound genera- 
tion. In fact, this is consistent with, Light hill's theory ' 


which shows that neither; energy- containing nor energy-dissipating 
eddies dominate the noise generation process; rather , noise gene- 


ration is dominated by an intermediate range of eddies: whose size 
is roughly eight times smaller than the energy- containing eddies. 
It has been shown elsewhere ^ that such intermediate 


scales of motion are already approximately Reynolds number inde- 
pendent at the Reynolds numbers of the numerical simulations. 

The principal source of Reynolds number dependence of laboratory 
jet flows is not the intrinsic viscous effect on the fluid motion, 
but rather is due to the viscous effect on the boundary layers 

at the jet exit which governs the flow through its initial and 

’ : . I 1 

boundary conditions. In the numer ica l simu latio ns , the viscous 
effect on ups h ream boundary and ini t i al cond i t ion s i s avo i ded by 
u se of mean profile s corr es pond i ng to Reyn o_l ds numbers of o rder 
80 0,000 not 20,000 . In this way it is expected that the numeri- 


cal simulations avoid most of the pitfalls of Reynolds number 
effects and thereby give results useful for the study of noise 
production. 

Some results of the numerical simulations are presented in 
Figs. 3-34. Here the three components of velocity are contoured 
within several slices of the computational domain at several 
evolution times (measured in seconds) and at. the three downstream 
locations D =4, 8, 12. These figures show the characteristic 

jet shape as well as indicating the nature of the motions re- 
solved by the simulations. 

The quadrupole moment distribution Q(x,t) for the throe 
simulations is contoured in Figs. 35-58, again for several times , 
locations, and slices. It is striking how the quadrupole moment 
distribution is so much more localized and spotty than the velo- 
city field itself, consistent with our treatment of a slab of 
the jet isolated from its surroundings. 


5. Discussion 


The results of the simulations discussed in Sec. 4 are being 
used together with Lighthill ’ s integral to determine the radiated 
sound field of the jet. For this purpose, Lighthill ' s integral 
is broken up into contributions from the several downstream slabs 

" f • • , . . . 

whose contributions are then superposed. The results of this cal- 
culation will be reported later. 

It is important to emphasize that noise generation .is possi- 
ble despite the periodic boundary conditions that are applied to 
the computational domain, because the cont ributions from neigh- 
boring slabs jars. distinct and do not cancel. 

There are at least three areas where the present computations 
are being extended and improved. They are: 

1) Inclusion of. a nonlinear eddy viscosity to model the 
effects of unresolved turbulence. These calculations with a sub- 
grid scale eddy viscosity are being compared with those with the 
molecular viscosity discussed here, in order to better gauge the 
Reynolds number effects on the simulations. 

2) Inclusion of compressibility, so that noise production 
at higher Mach numbers can be studied. 

3) Inclusion of realistic jet geometry, including jet : o 
growth, entrainment, and inflow-outflow boundary conditions . 
Because of limited spatial resolution on available computers, 
these flow simulations in realistic jet geometries are limited 
presently to about the first five jet diameters downstream. 

The latter improvement will permit the proper simulation of 
jets in the presence of a pressure gradient by allowing for jet 
expansion downstream. 

Another problem of interest jis the study of the present sim- 
ulations for evidence of coherent large scale flow structures, as 
observed recently by Crow & Champagne, Laufer, and others. Whi lc 
these large-scale flow structures probably do not dominate noise 
production, they may give a component of the noise that is con- 
trollable. In this way, it may be possible to achieve the design 
of minimum noise flows at a given flow rate and Mach number. 
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FIG. 3 - Contour plot of v 1 (x 1 ,x 2 ,x 3 > in the plane 

at X-4 diameters downstream from the jet exit. 

Here t=0, the contour minimum ' t ^ le contour 

maximum v max ~l»(2) and the contour interval Av=*l. (1) , 

The time t is the time of evolution with t=*0 the initial 
turbulent state. 





Fib. 4 - Same as Fig. 3 except the contours are for v 
in the plane X 2 =^-L. Here t=0, X«4, v ro ^ n =-l. 1 (2) , 
v -1.4(2), Av-l.tl). 



FIG. 5 - Same as Fig. 


in the plane x^^l. 


max 


! 9. (1) , Av=l. (l) . 


3 except the contours are for v 
Here t=0 , X=4 , v . =-9 . ( 1 ) , 

mm * ' » 






FIG. 7 - Same as Fig. 3 except the contours are for 
in the plane x^L. Here t=.020, X=4, v fflin »-l.l(2), 


max 


-1.1(2) , Av=l . (1) . 




FIG. 8 - Sane as Fig. 
in the plane x 2 =iL. 
v mAv=l-7(2), Av=l. (1) 


3 except the contours 

Here t=.020, X=4. v . . 

mm 


are 


for v 


- 1 . 2 ( 2 ) , 





FIG. 9 - Same as Fig . 3 except the contours are for v 

• i _ 

in the plane x_=~-L. Here t=.020, X=4, v . =-7 . ( 1 ) , 
z l min 

v_, =8 . (1) , Av=l . (1) . 






FIG. 10 - Same as Fig. 3 except the contours are for 
in the plane x 2 =^L. Here t=.020, X= 4 , v . =-6.4(1), 
v max =5 * 6(1) ' ^ v =8 . (0) . 




FIG. 11 - Same as Fig. 3 except the contours are for v 
in the plane x^-^-L. Here t=.031, X=4, v =-1.1(2), 
V n^ v =1 -4 (2) , Av=l. (1) . * mn 









FIG. 13 - Same as. Fig. 3 except the contours are for 

in the plane x,=k. Here t=.031, X=4, v' . =-9.(1), 

2 . 2 • min ' 

Av=l.(l). 




FIG. 14 - Same as Fig. 3 except the contours are for v 

in the plane x 0 =i-L. Here t=.031, X=4, v . =-8.(1). 

£ & mxn 

=8. (!) / Av=l. (1) . 








FIG. 15 - Same as Fig. 3 except the contours are for v 
in the plane x^jL. Here t=.051, X=4 , v . =-1.4(2) ^ 

v max =1 - 4 < 2 >' AV-1.(1). .min 
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FIG. 17 - Same as Fig. 3 except the contours are for ^ ' 

in the plane x 2 =i-L. Here t=.051, X=4, v . =-1.2(2). 

t t nun * • 

v mix =9. (1) , Av=l.(l). 


3 1 



FIG. 18 - Same as Fig. 3 except the contours are for v 
in the plane x -±L. Here t=.051, X*4, v . =-9 (1) 

min 

V m a «=l-K2). AV=1.(1). 















FIG. 20 - Same as Fig. 

in the plane x 2 =^-L. 

V =1.7(2) , Av=l . (1) . 
max 


3 excep 
Here t 


3 





*IG. 21 - 


« . y . j except the contours are 

V 2 in the plane x^L. | lere t . 0> jj, 8> y 

Av-1. (1) . min 



FIG. 22 - Same as Fig. 3 except the contours are for ^ 
in the plane x 2 = 2 L * Here t=0, X=8, v m i n = "8 . 1 (1) , 
=5.4(1) , Av=9 . (0) . 



FIG. 23 - Same as Fig. 3 except the contours are for 
V 1 in the plane Xj-jL. Here t=.020, X=8, v . =-1.(2 
V n^ v =1 -3(2) , Av-l.(l). 





FIG. 2 4 - Same as Fig. 3 except the contours are for'"*, 
V, in the plane x-=i-L. Here t=.020, X=8, v . *-1.1(2), 
v m „v = l-9<2), Av-l.(l). 






FIG. 25 - Same as Fig. 
v 2 in the plane x 2 =^-L. 
v m =1. (2) , Av=l . (1) . 


3 . 
















xcept the contours are for 
re t=0 , X=12 , v .=-1.2(2), 
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v max =6 • 3 (1) • ^v=9.(0). 
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*...v 




FIG. 31 - Same as Fig. 3 except the contours are 
for v x an the plane . Merc t«=.032, X=12, 

v min v max =1 * 6 (2), Av=i.(i). 
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FIG. 32 - Same as Fig. 3 except the contours are for^cT^ 
v ! in the plane x 2 =^L. Here t=.032, X=12, v mi *-l . l (2) 
V m =2. < 2 > » Av-1. (1) . 





33 • SdXHG Flrr *> 

v 2 in the plane x _ 1 L the c p n tours are £ or * 

* * a,ie X---L. Here t=.032 Y-1-) 

V -Q /i\ . , ^ X-12 , V . =-l 2/9% 

max 9 * Av«l.(i). min 
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FIG. 35 -Contour plot of the quadrupole moments in 
the plane x^^L at X=4 diameters downstream from the jet 
exit. Here t=0, the contour minimum is Q min =-1 . 4 ( 5) , 
the contour maximum is =4.(4) and the contour inter- 

JuclX 

val is AQ=1. (4) . 
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FIG. 36 - Sane as Fig. 35 except the contours are made i 

in the plane x_=-^L at X=4. Here t=0 , O . =-2.4(5), 
z z min 

Q,™=1.2(5) , AQ=2 . (4) . 




FIG. 37 - Same as Fig. 35 except the contours are made ^ 
in the plane x^^-L at X=4. Here t=0, Q m i n = "3 . 9 (5) # 
Q ina x=l-5 (5) , AQ=3. (4) . 




in i l ^ SXCept the contours are , 

in the plane x^L at 7-4. llere t=.020, Q 4 

Q max= 3 -6<5). AQ=4.<4). ml n *•* 


FIG. 39 - Same as Fig. 35 except 
in the plane x 2 =^'L at x= ^« here 
Q m , v =5 .6(5) , AQ=7 . (4) . 






FIG. 41 - Same as Fig. 35 except the contours are made * 
in the plane at X=4. Here t=.031, Q m i n = ~ 5 *( 5 )* 

Q„ =3. (5) , AQ=5 . (4) . 



fig. 42 - Sane as Fig. 35 except the contours are made 
in the plane x 2 =~l at X=4. Here t=.031, Q . =-1.5(6), 
AQ-1.(5>. . mn 




fig. Jj except the contours are made 
■L at X=4. Here t=.031 # o , =-1.2 (f>) 


nun 



FIG. 44 - Same as Fig. 35 except the contours are made"** 
in the plane x x =jL at X=4. Here t=.05i, Q min = _ l • 3 (6) » 
Q^ax' 8 * (5) * AQ-1. (5) . 








FIG. 46 - Same as Fig. 35 except the contours are made 

in the plane x,=-^L at X=4. Here t=.051, Q . 1 . 2 (6) , 

34 min 

Q max *3.3(6), AQ=3. (5) . 


FIG. 47 - Same as Fig. 35 except the contours are made 
in the plane at X=8. Here t=0, Q ra i n s “l» 6 (5) • 

Qjjjax* 8 * (4) , AQ=1. (4) . 






FIG. 48 - Same as Fig. 35 except the contours are made 

in the plane X=-L at X=8. Here t=0 , Q . =-3.6(5), 

^ 2 mm 

Q mav = l .5(5) , AQ=3.(4). 




FIG. 49 - Same as Fig. 35 except the contours are made 
in the plane x-pij-L at X=8. Here t=0, Q m ^ n =- 3 . 9 (5) f 
Qm*v=l-5(5) , AQ=3. (4) . 





FIG. 51 - Same as Fig 
in the plane x 2 =i-L at 
Q—-5.4(5>, AQ=9.(4) 


. 35 except the contours are maue 
X=8. Here t=.020, Q m - =-9.(5), 
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the contours are made**. 


Here t= 


>< <w- 1 - 4(5 >' 



FIG. 54 - Same as Fig. 35 except the contours are made 
in the plane x 2 =jL at X=12. Here t=0 , ^ m i n = “ 4 • 8 (S) # 

Q m =1.6 (5) , AQ=4 . (4) . 




FIG. 55 - Same as Fig. 35 ex 
in the plane at X=12. 

Qtn» v s l‘8(5) , AQ=6. (4) . 


the contours are made 


Here t* 


°' Q min=- 7 - 8 < 5 >' 













FIG. 56 - Same as Fig. 35 except the contours are made * 

in the plane x, =— L at X=12. Here t=.032, Q ■ =- 5.6(5) , 
i 2 min 

Qm av =4.9(5) , AQ=7 . (4) . 







FIG. 58 - Same as Fig. 35 except the contours are made 
in the plane X3=^-L at X=12. Here t=.032, n =”8. (5) , 
Qm*v B 4.8(5) , AQ=8. (4) . 




